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Abstract 

An approach is suggested for treating multiscale fluctuations in macromolecular 
systems. The emphasis is on the statistical properties of such fluctuations. The ap- 
proach is illustrated by a macromolecular system with mesoscopic fluctuations between 
the states of atomic orbitals. Strong-orbital and weak-orbital couplings fluctuationally 
arise, being multiscale in space and time. Statistical properties of the system are 
obtained by averaging over the multiscale fluctuations. The existence of such multi- 
scale fluctuations causes phase transitions between strong-coupling and weak-coupling 
states. These transitions are connected with structure and size transformations of 
macromolecules. An approach for treating density and size multiscale fluctuations by 
means of classical statistical mechanics is also advanced. 
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1 Introduction 



Macromolecular systems are assemblies of many atoms, ranging between hundreds or thou- 
sands to millions and many millions of atoms. Examples are viruses, ribosomes, and other 
assemblies, as discussed in Refs. 1-7 . Similar properties are also exhibited by artificial macro- 
molecules, such as nanoclusters 8-11 , quantum dots 12 ' 13 , and various ensembles of trapped 
atoms, whose description can be found in books 14-16 and review articles 17-29 . The nanosize 
objects can be employed in a variety of applications 14-29 , including such nontrivial ones as the 
creation of quantum dot biocomposite structures using genetically engineered bacteriophage 
viruses 30 . 

The systems of nanosizes possess peculiar properties by combining the features that are 
typical of macroscopic and microscopic systems. From one side, nanostructures contain 
many molecules and, hence, can be described by means of statistical mechanics, as applied 
to bulk matter. From another side, the finiteness of such structures can play an important 
role, making them dependent on their surrounding and influencing stronger fluctuations of 
observable quantities. 

One of the typical properties of complex macromolecular systems is the existence of 
processes involving widely separated time and length scales. Such systems exhibit two main 
types of fluctuations: fast fluctuations, related to atomic collisions and vibrations, and slow 
fluctuations, representing coherent motion of many atoms simultaneously. Between the short 
and long timescales, there is a gap that makes possible the distinction of these principally 
different kinds of motion. The slow coherent motion of many atoms is characterized by 
a broad range of timescales, which makes the description of such nontrivial dynamics a 
complicated problem. A multiscale analysis for describing the dynamics of nanosystems, 
involving Smoluchowski and Fokker-Planck equations, has been developed by Ortoleva et 
al. 1-7 . In the latter articles, one can find more information on the coherent multiscale 
motion of groups of correlated atoms. 

Phenomenological methods of treating multiscale mesoscopic fluctuations occurring in 
condensed-matter systems have been suggested by Frenkel 31 , Fisher 32 ' 33 , and Khait 34,35 . 

In the present paper, we also address the problem of describing complex macromolecu- 
lar systems displaying fluctuations in a wide range of space and time scales. However, our 
approach is different from that of Ortoleva et al. We aim at considering not the dynamics 
of multiscale fluctuations, but their statistics. Our approach is also different from the phe- 
nomenological methods of other authors 31-35 , as far as we aim at developing a microscopic 
theory. 

The realistic picture, we keep in mind, is as follows. Suppose, we observe the behavior 
of a system during the time interval that is larger than the typical fluctuation time of slow 
multiscale fluctuations, so that during this observation time, there arises a number of these 
fluctuations. This implies that the experimental measurements will produce the smoothed 
results averaged over all these fluctuations. Therefore, we do not need to know the detailed 
dynamics of the system at each moment of the experimental observation time, but we need 
to know only the smoothed picture averaged over all fluctuations that have happened during 
this time. This is what is called the statistics of fluctuations. In that sense, our approach is 
complimentary to that of Ortoleva et al. 1-7 . 

The paper is organized as follows. In Sec. 2, we delineate the general approach that will 
be employed in the following sections. In Sec. 3, the model of a macromolecular system 
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is derived, which we shall use for illustrating the approach. For generality, we consider 
a quantum model, keeping in mind that many biological systems to some extent display 
quantum properties 36-40 . In Sec. 4, the behavior of the system order parameters is analyzed. 
Of course, the use of quantum models is not always necessary. Therefore, in Sec. 5, we 
formulate the approach in the language of classical statistical mechanics. In Sec. 6, using the 
classical statistics, we show how multiscale density fluctuations can provoke sharp changes 
of a macromolecule size. Section 7 concludes. 



The aim of the present section is to give the ideological and mathematical foundation for 
the statistical approach of treating multiscale fluctuations. In Sec. 2.1, we describe how 
the characteristic space and time scales of statistical systems can be evaluated and explain 
in what sense multiscale fluctuations are mesoscopic. This picture makes it clear that the 
snap-shots of the system with mesoscopic fluctuations can be described by means of the 
manifold indicator functions, as is shown in Sec. 2.2. The specific feature of the mesoscopic 
fluctuations of being multiscale serves as a justification for the procedure of averaging over 
spatial configurations of such fluctuations, which is formulated in Sec. 2.3. This section also 
provides the general formulas that are necessary for considering the particular models of the 
following sections. 

2.1 Length and Time Scales 

Let us consider a macromolecular system consisting of many atoms or molecules. In order 
to give a precise classification of different types of motion and the related fluctuation types, 
let us, first of all, discuss the characteristic length and time scales typical of multiatomic 
systems. 

Atoms, or molecules, composing the system, interact with each other, which introduces 
the interaction length or scattering length, a s . The mean interatomic distance a is connected 
with the average atomic density p by the relation pa 3 = 1. Atoms can move without collisions 
at the distance of the mean free path A ~ 1/P a s- Groups of atoms are correlated and can 
move in a coherent way, if the group linear size is not smaller than the correlation length l c . 
The correlated atoms can coherently move in groups whose size is denoted as the fluctuation 
length If. The largest length is a linear size L of the whole system. The standard relations 
between these characteristic lengths satisfy the inequalities 



In many cases, it happens that a s is of order of A. The fluctuation length // is not just a 
single fixed quantity, but it can represent a dense set of values between l c and L. 

These characteristic lengths are connected with the corresponding characteristic times. 
The interaction time t int defines the duration of atomic interactions and is given by the ratio 
of the scattering length a s versus atomic velocity v ~ h/ma s , where m is atomic mass. The 
local equilibration time t\ oc is defined by the ratio of the mean free path A ~ 1/ pa 2 s over 
the atomic velocity v. And the correlation time t cor is given by the ratio of the correlation 



2 Microscopic Versus Mesoscopic Fluctuations 



a s < A < l c < If < L . 
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length l cor over velocity v. Summarizing, we have the interaction time 



a* ma 2 



tint — — — —jT- , (2) 



V 



local- equilibration time 



and the correlation time 



hoc — — — t > (3) 

v npa s 

tcor Z • (4) 

v a 

The time of experimental observation will be denoted by t exp . The typical relation between 
these times is 

tint *^ t\ oc < t cor < t eX p . (5) 

Though ti n t can be of order t\ oc . The coherent motion of atomic groups occurs during the 
fluctuation time tf that is between t cor and t exp , 

tcor ^ tj <C t eX p. 

Again, ty is not a single fixed quantity, but can densely fill the whole given temporal interval. 

Depending on the considered time interval, the system dynamics can be separated into 
several stages. The shortest is the interaction stage, corresponding to the time of atomic 
interactions, 

< t < t int (interaction stage) . (6) 

The second is the kinetic stage, when atoms have experienced just a few collisions, but have 
had yet not enough time to become well correlated, 

U n t < t < Uoc (kinetic stage) . (7) 

After the local-equilibration time, atoms become already correlated with each other and can 
form coherently moving groups of many atoms, 

Uoc <t< t cor (fluctuation stage) . (8) 

And, when the experimental observation time is sufficiently long, there exists the quasi- 
equilibrium stage, when a number of coherent fluctuations occur during the observation 
time, 

tcor < t < t exp (quasi — equilibrium stage) . (9) 

The stage is termed quasi-equilibrium, since, at each moment of time, the system is not at 
equilibrium, being nonuniform and dynamically non- stationary, but on the longer time scale, 
it can be treated as equilibrium on average, fluctuating around a quasi-stationary state. 

Estimating the above quantities for typical multiatomic systems of condensed matter 
type 1-7,19,41-43 , we have a s ~ a ~ A ~ 1CT 8 cm and l c > 1CT 7 cm. While for the typical 
times, we get t int ~ U oc ~ 10~ 14 s to 10~ 13 s, and t cor > 10~ 12 s. 

These estimates show that there exist two principally different types of motion. At atomic 
length and time scales, there are fast microscopic fluctuations related to atomic vibrations 
and collisions and characterized by the typical times t int ~ 1CT 14 s. And there are slow 
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fluctuations corresponding to the coherent motion of correlated atomic groups, with the 
typical times t cor > 10 -12 s. The slow fluctuations occur in the diapason of lengths, If, and 
times, tf, satisfying the inequalities 

A < l c < If < L , t loc < t cor < tf < t exp . (10) 

In that sense, these coherent fluctuations can be called mesoscopic. 

As has been stressed by Ortoleva et al. 1-7 , there are two crucial points distinguishing 
these two types of fluctuations. First, microscopic fluctuations of individual atoms and 
mesoscopic fluctuations of correlated atomic groups are separated by a large gap. Thus, the 
former are characterized by the time scale of 10~ 14 s, while the lowest scale of the latter is of 
order 1CT 12 s. Second, the slow mesoscopic fluctuations do not have just one time scale, but 
they can happen in a wide temporal range, starting from 1CT 12 s and up to the observation 
time. That is, the mesoscopic fluctuations exhibit multiscale motion. 

2.2 Snap-Shot of Heterogeneous State 

Let us assume that we can make a snap-shot of the system at a fixed moment of time. In 
the presence of mesoscopic fluctuations, the system is non-uniform, being separated into 
several regions with different properties. These properties can be characterized by order 
parameters or order indices 44,45 . Actually, the mesoscopic collective fluctuations are nothing 
but the slow fluctuations of the order parameters 1-7 . Thus, at each moment of time, the 
system can be highly nonuniform, which is termed heterogeneous state. Generally, the order 
parameters can be associated with some phases, because of which this heterogeneous state, 
with nonuniform order parameters, can also be called heterophase state 41-43 and the system 
parts can be termed the phases. 

Making a snap-shot of a heterogeneous system, one can distinguish the spatial regions 
with different order parameters. Then the whole system volume can be represented as the 
union of these spatial manifolds: 

V V 

with the related manifold measures 

V = mesV , V v = mesV„ , (12) 

where the index v = 1,2, .. . enumerates the phases with different order parameters. The 
separation into the subregions can be done by invoking the notion of the equimolecular 
separating surface, when the total number of atoms N is the sum of atoms in each of the 
subregions, 

N = Y,N V . (13) 

V 

The definition and properties of the equimolecular separating surface have been introduced 
and described by Gibbs 46 ' 47 and are widely employed in chemical and physical literature 48 ' 49 . 
This separating surface is defined in such a way that the surface density be zero, because of 
which the additivity of the atomic numbers is exact. Another convenience of employing this 
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separating surface is that the volume and all extensive thermodynamic characteristics of the 
whole system are additive with respect to the volumes and thermodynamic characteristics 
of its parts 46-49 , similarly to property (13). The equimolecular separating surface can be 
defined for any spatial region of arbitrary size 46-49 . 

Mathematically, the separated spatial regions can be described by the manifold indicator 
functions 50 . If the phase, labeled by the index v, occupies the volume V,, the related 
manifold indicator function is defined as 

«'>-{£ (14) 

The set 

£ = {^(1-): i/ = l,2,...; r G V} (15) 

of all indicator functions uniquely describes the given spatial configuration. 

The system with the fixed spatial configuration of local phases should be described by the 
Gibbs local-equilibrium ensemble 46 ' 47 . For this purpose, it is necessary to define the local 
single-particle energy operators K v (r) and the local interaction energy operators $(r,r'). 
The Hamiltonian of the //-phase is written in the form 

= J [k u (v] + E„(r)] £„(r) dr + l -J ® v {r,r')Ur)U*') dvdr' , (16) 

in which the integration is over the whole system volume (11) and £„(r) is an operator of 
external fields. If there are no external fields, the term £„(r) plays the role of a source mim- 
icking the local influence of the separating surface. Characterizing just local perturbations 
near the separating surface, this term has negligible influence, as compared to the single- 
particle bulk term K v (r) and it can be set to zero after calculating observable quantities, as 
it is done in the Bogolubov method of quasi-averages 51 . However, it can be essential if there 
are external fields or the system volume is small. 
The total system Hamiltonian 

#(O = 0#,(O (17) 

V 

is the direct sum of terms (16). Each of the terms H v (£) acts on a weighted Hilbert space 
1-L U and the total Hamiltonian (17) is defined on the fiber space 

H = ®H V , (18) 

V 

which is the tensor product of the weighted Hilbert spaces 43 . According to the definition 
of the equimolecular separation 46,47 , preserving the additivity of extensive quantities, the 
operators from the algebra of local observables are represented as the sums 

MO = (&MO, (19) 

V 

with A v {£) acting on % v and on space (18). 
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In the present section, for generality, we describe the approach that is valid for quantum 
systems. As has been mentioned in the Introduction, many biological systems, even of rather 
large sizes, exhibit quantum properties 36 " 40 . Then, in Sees. 3 and 4, we consider a quantum 
model illustrating the approach. In Sees. 5 and 6, we explain how the approach can be 
used for characterizing the size transformations induced by multiscale density fluctuations 
in classical systems. 



2.3 Averaging over Multiscale Fluctuations 

Any experiment, accomplished over the time period essentially larger than the correlation 
time (4), will produce the results corresponding to the smoothed picture averaged over many 
mesoscopic fluctuations that occurred during the observation time. This implies that we do 
not need to follow the complicated dynamics of all those fluctuations that have happened 
during the quasi-equilibrium stage (9), but we have to understand how to describe the related 
averaged picture. Multiscale fluctuations are characterized by a wide diapason of time and 
space scales 1-7 . Therefore averaging over them requires to consider various possible snap- 
shot configurations described in the previous subsection. In mathematical terms, this means 
that it is necessary to define functional integration over the manifold indicator functions (14). 
Such a functional integration has been rigorously defined in the previous papers 41-43 ' 52 ' 53 , 
where all mathematical details can be found. Not going into these details, we denote the 
differential measure of the functional integration over sets (15) of manifold indicator functions 
(14) as Vi. 

First of all, we have to understand how the functional integration over sets (15) can be 
realized for the polynomial functionals of the form 

C„(0 = Yl Y ■ ■ ■ Y J C w-<>n( T u r 2, • • • , r„)^ 1 (r 1 )^ 2 (r 2 ) . . . f„„(r n ) dr l( lr 2 ...dr n , 



V\ V2 v. 



(20) 

where the integration is over the system volume (11). Here the kernel C VlV2 _ Mn (•) is an 
arbitrary operator function. When functional (20) represents observable quantities, the 
kernel function is self-adjoint. 

The following proposition 42 ' 43,52,53 is valid. 

Theorem 1. The functional integration of the polynomial functional (20) over the 
manifold indicator functions (14) gives 



C n (0 Vi = C n (w) , (21) 
where 

C n {w) = YY---Y W ^ W ^--- W ^ J C ^2...u n (r 1 ,r 2 ,. . . ,r n ) dr 1 dr 2 . . . dr n , (22) 



and 



V 

is the geometric weight of the v— phase. □ 



w v = \-\ £,(r) dr (23) 
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The operators of observable quantities act on the Hilbert fiber space (18). Each of the 
fiber sections H v contains information on the degrees of freedom corresponding to microscopic 
atomic fluctuations. The variables corresponding to mesoscopic multiscale fluctuations are 
represented by the manifold indicator functions (14). A macromolecular system, with these 
two types of degrees of freedom is characterized by the Gibbs quasi-equilibrium ensemble 46,47 . 
It is possible to introduce temperature T as a thermodynamic parameter corresponding to 
the temperature of thermostat surrounding the given system. In a more general interpreta- 
tion, temperature can be associated with the intensity of noise, caused by random external 
perturbations 29 . The related thermodynamic potential reads as 



J Tc H e- fiH ^ Vi , (24) 



F — —T In J Tr n e 

where the trace is taken over the microscopic degrees of freedom and (3 = 1/T. Here and 
everywhere below, we use the system of units where the Planck constant and Boltzmann 
constant are set to one (h = 1, ks — 1). The following statement holds 42,43,52,53 . 

Theorem 2. The integration over the manifold indicator functions (14) in the thermo- 
dynamic potential (24) gives 

F = -ThiTr u e-t S = Y,F V \ ( 25 ) 

V 

here the renormalized effective Hamiltonian is 

i? = 0if„, (26) 

in which 

o 

H v = w u (K v + £„J + $„ , 

K v = J K u (r) dr, E„ = J t u (r) dr , $ v = J $,(r, r') drdr' , (27) 
the partial thermodynamic potentials are 

F^-TlnTr^e-^ , (28) 
and the geometric weights w u are the minimizers of the thermodynamic potential (25), 

F — abs mm F({w v }) , (29) 

{w v } 

under the conditions 

^uv = l, 0<w v <l, (30) 

V 

defining the set {w u } of the geometric weights as a probability measure. □ 

The observable quantities are given by the expectation values of operators (19) from the 
algebra of local observables. These expectation values 

(A) = J TruKOMO n (31) 
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contain the trace over the microscopic degrees of freedom and the averaging over multiscale 
mesoscopic fluctuations with the statistical operator 

m _ exp{-/3tf (O 



fTr n exp{-PH(Z)}Vf 

It is possible to prove the following theorem 42 ' 43 ' 52 ' 53 . 

Theorem 3. Integration over the manifold indicator functions in the expectation value 
(31), for large N >> 1, yields 

<i> = 5>^>' (32) 



with the terms 
partial statistical operators 

and the notation 



(A v ) = Tr H JvA v , (33) 

_ exp(-PH u ) 
Tr Hv exp(-pH v ) 

A u = lim A v {£) (35) 

implying that all £„ are replaced by w„. □ 

Sometimes, it is more convenient to resort to the equivalent notation 

(A v ) = (A) u , (36) 

which we shall use when appropriate. 

Among the observable quantities, there are the order operators fj u , whose averages define 
the order parameters 

Vu = (fj v ) (37) 

allowing us to distinguish different phases. 

As a particular example, let us consider the case of two phases, when v — 1,2. Then it 
is convenient to invoke the notation 

Wi = w , w 2 = 1 — w , (38) 

explicitly taking account of normalization (30). The minimization condition for the thermo- 
dynamic potential (25) becomes 

dF d 2 F 

dw ^ ' dw 2 > ^ ^ 
The first of the above equations results in 
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While the inequality in Eq. (39) becomes 

d 2 H 




(41) 



In view of the form of Hamiltonian (26), Eq. (40) reduces to 



$ 2 + K 2 - 

W = »,+*, • (42) 

where 

K V = K V + V V , K v = (K v ) , 

E„ = (£„) , = ($„) . (43) 

Equation (41) has the meaning of the stability condition and reads as 

•' + *> > "((^))- (44) 

Since the right-hand side here is positive, the necessary stability condition is 

$i + $ 2 > . (45) 

This tells us that multiscale mesoscopic fluctuations arise not in an arbitrary system, 
but only in those systems where the stability conditions are satisfied. The stability con- 
ditions impose restrictions on the system parameters as well as on the external param- 
eters, defining the regions of admissible values of these parameters where the multiscale 
fluctuations can exist. The concrete constraints, following from inequalities (44) and (45) 
essentially depend on the considered model. The possibility of occurrence of mesoscopic 
fluctuations has been analyzed, for several models of condensed matter, for instance, for an- 
harmonic crystals 52 , lattice-gas model 54 , Hubbard model 55 , Vonsovsky-Ziener model 56 , and 
high-temperature superconductors 57-59 . Somewhat similar situation happens for the model 
of mixed nuclear matter 60 ' 61 . 



3 Model Macromolecular System 

In the present section, we shall derive a model of a macromolecular system of N atoms, with 
iV assumed to be large. This model will be used for illustrating how the approach of Sec. 2 
works. For generality, we consider a quantum system, since many biological systems to some 
extent display quantum properties 36-40 . Typically, quantum effects in biomolecular systems 
come from fluctuations of protons and long-range electron transport /photo excitation in 
macromolecules. Of course, quantum effects not always are important. And in Sections 5 
and 6, we show how the approach could be used for describing multiscale fluctuations in 
classical systems. 

Meanwhile, for generality, we consider a generic quantum system described by field op- 
erators ip(r). The standard expression for the single-particle energy operator is 

K = I V> f (r)#(r)V>(r) dr , (46) 
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with the single-atom Hamiltonian 

H(r) = - ^ + U(r) , (47) 

where U (r) is an effective potential making the system confined. Additional energy, caused 
by external fields, corresponds to the operator 

t = J ^t(r)/i(r)^(r) dr . (48) 

And the interaction term has the standard form 

<f> = J ^(r)^(r')<S>(r - r')^(r / )^(r) drdr' , (49) 

with $(r) being the pair- interaction potential. 

For what follows, we need a basis of orthonormalized wave functions, which can be chosen 
as the set of the eigenfunctions of the Schrodinger equation 

H(r)ip n (r - Tj) = E nj ip n (r - rj) , (50) 

defining the localized orbitals 44 centered at rj, with j = 1, 2, . . . , N. Being orthonormalized, 
these functions enjoy the property 



J V4(r - ri)V>n(r - Vj) dr = 



Because the considered system is confined, its spectrum is discrete. 

The field operators can be expanded over the basis of the localized orbitals: 



V>(r) = 5^c ni ^„(r-r j ) . (51) 



nj 



Substituting this expansion into Eqs. (46) yields the single-atom term 

K = Y,E n3 cifn, ■ (52) 

nj 

And using the expansion in Eq. (48) results in 

t = E E h T c liC nj , (53) 

mn ij 

where 

= j PJt - rOMr)^n(r - r,-) dr . (54) 
Respectively, for the interaction term (49), we get 

$ = V" $«l«2n3Tl4 t t „ . „ . / 55 \ 

^ J1J2J3J4 C riUl C ri 2 J2 Cn 3J3 C ri 4 J4 • 1< J '-V 

nynwzm jxjijzjA 
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In order to take into account that at each spatial location rj there is just one atom, we 
impose the unipolarity constraint 

^ ] C nj C nj = 1 j CmjCnj = . (56) 



n 



And the condition that atoms are localized reads as 

C rni C nj = 3ij C mj C nj • (57) 

Conditions (56) and (57) make it straightforward to simplify the above expressions, so 
that Eq. (53) reduces to 

^ = EE^r4^ (ss) 

ran j 



and Eq. (55) becomes 



* = E E y ~' n ' c L4 c ^ c "'- (59) 

nin2n3ri4 i^j 



where the notation 



jrmnrnn _ ^mnra n _i_ ^r. 



is used, with the sign plus or minus, depending on Bose or Fermi statistics of atoms, respec- 
tively. 

Assume that only two lowest energy levels, defined by Eq. (50), are occupied, so that 
the index enumerating the levels is n — 1, 2. The wave functions, corresponding to different 
energy levels, usually possess different symmetry, which makes zero some of the matrix 
elements y^ nm n . The remaining matrix elements can be combined into the expressions 

A, ee \ K 111 + V™ 22 + 2^ 221 ) , B i3 = \ ( V j™ + - 2V^) , 
Also, let us introduce the notations 

3> = 2^ £(^y + 3y) , (61) 

ee £ 2j + hf 3 - - fc}} + E C« > ( 62 ) 

and 

= - Jig . (63) 

It is convenient to accomplish an operator transformation by introducing the pseudospin 
operators 



(64) 
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These operators enjoy the spin algebra independently from whether the considered atoms 
are bosons or fermions. The inverse transformations are 

c \. C2j S-j iS? , 4 jClj = S* + iS? . (65) 

By the above definition, the operator S* describes the state imbalance between the or- 
bitals, the operator Sj corresponds to the interstate current, and characterizes transitions 
between the orbitals. With these notations, we have 

K + ± + Q = NE -Y f {VjS* + BjS*) + £ Q A v + B ^ S J ~ WS£ 

3 i+3 

Note that operators (64) are called pseudospin, since, though they enjoy the spin al- 
gebra, they do not represent real spins, but just serve as a convenient mathematical tool 
characterizing atomic transitions between orbitals. 



4 Statistics of Multiscale Fluctuations 

To take into account the influence of multiscale fluctuations, we follow the approach presented 
in Sec. 2. The main order parameter distinguishing different phases that exhibit multiscale 
fluctuations is the average strength of interorbital transitions 

*^£W>"- ( 66 ) 

3 

This quantity characterizes the intensity of interorbital transitions, or the strength of in- 
terorbital coupling. 

Mesoscopic multiscale fluctuations are known to be connected with the molecular struc- 
ture, molecular shape and size 1-7 . These features are incorporated into the value of the 
order parameter (66). If there are no transitions between the orbitals, then the system is in 
its ground state corresponding to the first orbital with the lowest energy. In the majority of 
applications, the ground-state orbital can be well approximated by a Gaussian. But when 
there appears the coupling between the orbitals, the higher orbitals starts contributing to 
the values of all observable quantities. The higher orbitals possess the larger mean radius. 
Hence, the involvement of these orbitals implies the increase of the system size, such as 
swelling. Therefore a sharp variation of the order parameter (66) is related to an abrupt 
swelling (or squeezing) of the macromolecule. Generally, one can consider a variety of shapes 
and sizes. Here, for simplicity, we take into account two possible phases enumerated with 
v — 1,2, defining their relation according to the inequality 

si > s 2 (67) 

between their order parameters. 

Here we are developing a general approach for describing statistics of multiscale fluctua- 
tions. And we aim at demonstrating that multiscale fluctuations can lead to the occurrence 
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of sharp phase transitions. We treat a model system that is convenient for such an illus- 
tration. It is not our aim to study particular substances, such as viruses. Their detailed 
investigation cannot be done in the frame of one paper, but should be a topic for separate 
publications. 

In the case of the macromolecular model of Sec. 3, following the method described in 
Sec. 2, we get the effective renormalized Hamiltonian 

H = H 1 @H 2 , (68) 

which takes into account the existence of mesoscopic multiscale fluctuations, and where 

H v = w u NE -w„J2 (^ + Btf) + £ (5 Ai i + B ^ SX 3 - ^ S i) ■ ( 69 ) 

According to expressions (60), the value of By is much smaller than that of iy, hence, can be 
neglected, to a first approximation. Also, keeping in mind a large molecule of many atoms, 
it is admissible to consider the parameters Qj and Bj as weakly depending on the atomic 
number, replacing them by Q and Bq, respectively. 

Accomplishing calculations with Hamiltonian (68), one meets the standard problem of 
the necessity to decouple atomic correlation functions. The simplest kind of decoupling 
is the mean- field approximation that, however, does not take into account the interatomic 
correlations, which then are completely lost. A much more elaborated is the Ter Haar 
approximation 62 that preserves in full pair correlations. The disadvantage of the latter ap- 
proximation is its high complexity. The intermediate situation is provided by the Kirkwood 63 
approximation that is sufficiently simple, at the same time taking account of atomic corre- 
lations. Employing the Kirkwood approximation in our case, we have 

(S?Sl) v = <%(S?) v (Sl) v , (70) 

where g v - is a pair correlation function for the //-phase. In terms of operators, acting on the 
weighted Hilbert space % V) this is equivalent to the transformation 

S?S? = g\ 3 ((S?)„Sf + S?(Sf) v - (S?)„(S;) V ) . (71) 

As is evident, averaging Eq. (71) yields exactly Eq. (70). 

Implementing the above approximation, we will need the notations 

Also, it is convenient to define the correlation parameters 

Then, Hamiltonian (69), in the Kirkwood approximation, takes the form 

(2 2 \ 

w v E + ^A + ^- Jg v slj - K fi5 7 + KA) + wlJg v s v )S?] . (74) 
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In order to simplify the formulas, we introduce the dimensionless reduced Gibbs potentials 

B 4' G ^m < 75 > 

and the dimensionless quantities 

u = — , w= j' J ' ' ' 

Also, in the following expressions, temperature will be measured in units of J. Then for the 
system Gibbs potential, we get 

G = G 1 + G' 2 , G„ = e„- Tin (2 cosh , (77) 

where 

e u = w u -j + -j- (2m + ^s^) , fij, = a/w 2 + (/i + w u g u s u ) 2 . (78) 
Except the order parameters (66), we can calculate the mean orbital imbalance 

The mean interorbital current is found to be zero, 

(Sj) = , (80) 
as it should be for a quasi-equilibrium system. The orbital imbalance (79) reads as 

(81) 

And for the main order parameters, the orbital coupling (66), we find 

h + w u g u s v (wA\ 

s »= o, v Unh {-zr)- (82) 

In view of expression (62), the tunneling parameter u, defined in Eq. (76), is small, w < 1, 
because of which the orbital imbalance (81) is always small, x u <^1. This confirms that the 
main order parameter is the mean orbital coupling (66), for which we have expression (82). 
In the Gibbs potential (77), the sum 

w u E = E = const 

V 

reduces to a constant, hence can be omitted. Keeping in mind that w C 1, we obtain 



,2 



•'«":; . . . , / w u h + w v g v s l 



G v = -f-(2u + g„sl)-Tln 



2 cosh 



2T 



(83) 
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And the orbital-coupling order parameter (82) becomes 



w = L ' ; ~J ~ 2 • (85) 



^^( ■^ + »' ) . (84) 

As in Eq. (38), we set w\ = w. Then, minimizing the Gibbs potential (77) with respect 
to w, we obtain 

2u + h(si - s 2 ) - g 2 s\ 
Au - gis\ - g 2 s\ 

which is the representation of form (42) for the considered model. Expression (85) is, actually, 
the additional order parameter defining the geometric weight of the phase with strong orbital 
coupling. 

Thus, we have the system of interconnected equations for the order parameters si, s 2 , and 
w. The phase with the strong orbital coupling is connected with essential atomic correlations, 
because of which the correlation parameter (73) for this phase can be set as g\ = 1. While 
the phase with the weak orbital coupling is characterized by the correlations only between 
the zo nearest neighbors, when the correlation parameter (73) is 

g 2 = ^ ~ (N » 1) . 

We solve the system of equations for s±,s 2 , and w numerically for different parameters 
u and h. The parameter u — A/J describes the ratio of repulsive to attractive interactions, 
because of which it can be called the repulsion parameter. And the parameter h measures 
the strength of an external field acting on the system. When, solving the equations, we get 
several solutions, then that of them is to be chosen, which provides the absolute minimum 
for the Gibbs potential (77), that is, which corresponds to the most stable system. Also, 
constraint (67), distinguishing the phases, is to be valid. The latter always holds for all T if 
h — 0, when s 2 = 0. If h is finite, there can arise a small region of T <C h in the vicinity of 
T = 0, where s 2 is not well defined. In that case, we define it by an analytical continuation 
from the case of zero to finite h, by employing the linear extrapolation typical of numerical 
analysis 64 . Strictly speaking, mesoscopic fluctuations might be absent. Therefore, looking 
for the most stable solution, we have to compare the Gibbs potential (77), characterizing a 
macromolecule with multiscale fluctuations, and the Gibbs potential G(w = 1) corresponding 
to a macromolecule without mesoscopic fluctuations, when w = 1 and there is just one order 
parameter given by the equation 



Si = tanh 



h+_si 
2T 



The most stable solution corresponds to the minimal of all those Gibbs potentials. The 
results of numerical calculations are presented in Figs. 1 to 5. 

Depending on the value of the parameter u, it is possible to separate several regions with 
qualitatively different behavior. First of all, when u = 0, mesoscopic fluctuations do not 
arise, as well as they are suppressed for very strong external field h — > oo. Therefore, we 
shall concentrate our attention on the region of u > and not too large external field h. We 
analyze the behavior of the order parameters w, si, and s 2 as functions of T for varying h. 
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In the region < u < 0.5, the occurrence of multiscale fluctuations leads to the appear- 
ance of first-order phase transitions. Figure 1 shows the order parameters w, si, and S2 as 
functions of temperature T for different values of the field h, at fixed u = 0.3. At zero T, 
there are no mesoscopic fluctuations, but they arise for T > 0. 

Larger values of the repulsion parameter u make it possible, for some fields h, the ap- 
pearance of mesoscopic fluctuations even at T = 0. This is illustrated in Figs. 2 and 3. 
The difference of Fig. 3 from Fig. 2 is in the existence of a region of external fields, where, 
instead of a first-order phase transition, there occurs a sharp crossover. 

When u > 1.5, sharp crossovers are presented for low T, while the first-order transitions 
appear at higher T. This is shown in Fig. 4. 

The first-order phase transition lines are presented in Fig. 5, as a function of h, at fixed u, 
and as a function of u, at given h. Increasing the external field h shifts the phase transition 
point to the right, while increasing u diminishes the transition temperature. When u — > 0, 
then the transition temperature tends to infinity, in agreement with the fact that at u — 0, 
there are no mesoscopic fluctuations. For example, at h — 0.5 and u = 0.001, the transition 
temperature is T = 626. 

The general understanding is as follows. Taking into account the existence of mesoscopic 
multiscale fluctuations is crucially important. Their existence leads to the occurrence of first- 
order phase transitions or sharp crossovers between the strong-orbital coupling and weak- 
orbital coupling states of macromolecules. Increasing the intensity of noise, or temperature, 
results in the appearance by a jump of mesoscopic fluctuations, when the geometric weight 
of the strong-coupling phase, w, abruptly falls down. Respectively, the strong-coupling 
order parameter si also falls down by a jump, while the weak-coupling order parameter $2 
suddenly increases. The point of such a phase transition depends on the values of the system 
parameters u and h. 



As has been mentioned above, we have considered the case of multiscale fluctuations in a 
quantum macromolecular system. This has been done for generality, keeping in mind that 
many biological systems possess quantum properties. However, multiscale fluctuations can 
equally arise in classical systems. In order to show how they should be described in that 
case, we present below the approach for considering the statistics of multiscale fluctuations 
in the language of classical statistical mechanics. 

Let us consider a system composed of N atoms or molecules, enumerated with the index 
i = 1, 2, . . . , N. The collections of all spatial variables and momenta are denoted through 
the sets 



5 Classical Multiscale Fluctuations 




P N = {Pi,P2,---,Pjv} , 



(86) 



N 



N 



dr N = l[dr 



dp N = l[dp t . 



i=i 



i=l 
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Kinetic energy of iV atoms is 

N 2 
i=l 

And the potential energy writes as 

1 1 N 

2$n = 2E^- r ^ ( 88 ) 

The overall consideration is the same as in Sec. 2. Under a fixed spatial separation of 
phases, we have the Hamiltonian H(£) = J2 U H U , with the terms 

N 2 -i N 

1=1 ~ " i^jr' 

where £ vi = ^(r,) is the manifold indicator defined in Eq. (14). The latter is normalized as 

N 



J2U^) = K, (90) 



i=i 

with jV„ being the number of atoms composing the z/-th phase. Accomplishing the averaging 
over the configurations of multiscale fluctuations, according to Sec. 2, we get the effective 
Hamiltonian, similar to Eq. (26), with the terms 

H v (r N ,p N ) = w u K(p N ) + ^ $(r*) , (91) 
where the phase weights are 

= ^ ■ (92) 

The canonical partition function for a heterophase system with multiscale fluctuations 
reads as 

dr N dp N 



Z(T, V,N) = H Z U (T, V, N) , Z„(T, V,N) = J exp {-f3H u (r N ,p N )} m(2]r)3N , 

(93) 

where the integration is over the phase volume 

Q = V N x M. 3N = {i-j G V, pi G R 3 : i — 1,2, ... , N} . (94) 
Integrating out the momenta, we have 

^i-(S w H-^1f (95) 

Having defined the system partition function, it is straightforward to calculate the desired 
thermodynamic characteristics. 
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6 Density and Size Fluctuations 



To be more specific, let us consider the case when multiscale fluctuations correspond to 
mesoscopic density fluctuations, so that, at each moment of time, the system consists of 
randomly distributed mesoscopic regions having density pi, which are intermixed with the 
regions of density p 2 . And let, for concreteness, the first phase be more dense than the 
second one: 

Pi > p 2 • (96) 

The densities p v play here the role of the order parameters distinguishing different phases, 
the dense phase (u — 1) and the rarified phase {y = 2). 

It is clear that, when the majority of N atoms are in the dense phase, the system volume 
V\ ~ N/pi is smaller than the volume V2 ~ N/p2, when almost all atoms are in the rarified 
phase. The transition from the smaller volume V\ to the larger volume V2 implies the system 
swelling, while the opposite transition from the larger volume V2 to the smaller volume V\ 
describes the system squeezing. 

Since in this consideration, the system volume is not conserved, it is necessary to employ 
the isobaric ensemble, for which the system partition function is a function of temperature 
T, pressure P, and the total number of atoms N. For the heterophase system, the isobaric 
partition function is 

Q(T, P, N) = TT Q U (T, P, N) , Q U (T, P,N)= / e~ ppv Z„(T, V, N) dV . (97) 

So that the Gibbs potential is the sum 

G(T,P,N) = J2G„(T,P,N) (98) 

V 

of the terms 

G U (T,P,N) = —Tin Q U (T, P, N) . (99) 

The phase weights are defined by minimizing the Gibbs potential over w u , under the 
normalization condition w\ + W2 = 1. Setting w = wi, W2 — 1 — w makes it possible to write 
the minimization condition as 

BC 

_ = , G = G(T,P,N). (100) 

The observable quantities can be represented as the averages over the distribution func- 
tion 

U(r N ,p N ,V) = i- exp{-0(H v + PV)} , 

v 

where the Hamiltonian H u = H u (r N ,p N ) is given in Eq. (91). This distribution is normalized 
by the condition 

h \lo ' P • 'WW 

The average of a function A u = A u (r N ,p N ) is defined as 



dV=l . 
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For example, the mean kinetic energy of the z/-th phase, takes the form 



which can be reduced to 

3TN 
2w v 

And the mean potential energy of the v-th phase is (1/2)$^, where 

Minimizing the Gibbs potential (98), according to condition (100), with taking into 
account that 

dG I dH v ' 



we find the equation for the weight of the dense phase 

$ 2 + K 2 - K x 



w = 



$1 + $2 



which is similar to Eq. (42). 

The total system volume is defined through the derivative 

dG 

V(T,P,N) = — = V 1 + V 2 , 

with the partial volumes 

f)G 

K = ^ = K(T,P,N). 
In this way, we find the densities of the competing phases 



N v ( dG„\ 1 

Pv = — = w v N 



At low temperature, when w\ > w 2 , almost all atoms are in the dense phase, so that iV ~ 
Ni and V ~ Vi, though the occurrence of the multiscale density fluctuations results in the 
existence of the regions of the rarified phase. Rising temperature leads to the intensification 
of the mesoscopic density fluctuations and to the increase of the weight w 2 . Then, these 
multiscale density fluctuations provoke the transition to the state, where w 2 > wi, so that 
the system volume increases, becoming close to V 2 , though a small admixture of the dense- 
phase fluctuations can remain. Such a swelling transition is illustrated by Fig. 6. 

7 Conclusion 

Macromolecular systems are investigated experiencing, in addition to fast atomic fluctua- 
tions, slow multiscale fluctuations corresponding to the coherent motion of groups of many 
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atoms. In the process of this motion, the system symmetry can be locally broken or restored 
in a spontaneous way 65 . When the experimental observation time of the system is sufficiently 
long, being much longer than the atomic correlation time, one does not need to study the 
detailed dynamics of all atoms. In such experiments, one observes the smoothed picture 
averaged over many mesoscopic multiscale fluctuations. Then what one has to know is the 
statistics of multiscale fluctuations, that is, the average influence of these fluctuations on the 
values of observable quantities. 

A method is suggested for treating the statistical influence of multiscale fluctuations on 
observables. The method is based on the Gibbs theory of quasi-equilibrium systems, with 
averaging over the ensemble of multiscale fluctuations. As a result, one comes to an effective 
statistical system with renormalized interactions and with the influence of the multiscale 
fluctuations incorporated into geometric weights of competing phases. 

The method is illustrated for a macromolecular system exhibiting two possible states 
differing by the strength of coupling between atomic orbitals. Multiscale fluctuations in 
space and time occur between these two states, of weak-orbital and strong-orbital coupling. 
The order parameters of the system are calculated. The presence of mesoscopic fluctuations 
is shown to lead to sharp phase transitions between the phases of orbital weak-coupling 
and strong-coupling. Such phase transitions can be accompanied by an essential swelling or 
squeezing of macromolecules. 

As is shown, the existence of mesoscopic multiscale fluctuations may result in drastic 
changes in the properties of macromolecules. When varying temperature, or other thermo- 
dynamic parameters, macromolecules can change their properties not gradually, but in a 
sharp phase transition. This is necessary to take into account in analyzing the behavior of 
macromolecules under the variation of external conditions. For example, considering the 
problem of the origin of life on Earth, one assumes that life appeared after the primordial 
ocean was cooled down 66 ^ 70 . In this respect, one often argues that it looks strange that 
biomolecules have arisen so suddenly, while the process of cooling was gradual. However, 
this does not seem strange if we take into account mesoscopic multiscale fluctuations. As 
is shown above, under a slow gradual cooling, the properties of macromolecules can change 
by an abrupt jump. It is, probably, the occurrence of mesocopic fluctuations that facilitated 
the spontaneous generation of life on Earth. 

The aim of the present paper has been threefold. First, to develop a general approach 
allowing one to describe the statistics of multiscale fluctuations in macromolecules, without 
the need of studying their dynamics. Second, to demonstrate the use of the approach by 
a model that, though being simple, nevertheless exhibits rather rich properties, with phase 
transitions provoked by the multiscale fluctuations. Third, to emphasize that the approach 
is equally applicable to quantum as well as to classical systems. We have limited ourselves 
by these principal problems. Applications to particular biomolecules, having complicated 
structure, require separate investigations, with extensive computer calculations that are out 
of the scope of the present paper and are planned for future research. 
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Figure Captions 



Fig. 1. Macromolecule order parameters: (a) geometric weight of the strong-coupling 
phase w; (b) strong-coupling phase order parameter s±] (c) weak-coupling phase order pa- 
rameter s 2 . Dependence on dimensionless temperature is shown for the repulsion parameter 
u = 0.3 and varying strength of external field: (1) h = 0.01; (2) h = 0.1; (3) h = 0.2; 
(4) h = 0.3; (5) h = 0.5; (6) h = 1. The corresponding first-order transition temperatures 
are: (1) T = 0.145; (2) T = 0.241; (3) T = 0.336; (4) T = 0.436; (5) T = 0.660; (6) 
T = 1.422. 

Fig. 2. Macromolecule order parameters: (a) w; (b) Si, (c) s 2 . Dependence on dimen- 
sionless temperature is shown for the repulsion parameter u = 0.6 and varying strength of 
external field: (1) h = 0.01; (2) h = 0.1; (3) h = 0.2; (4) h = 0.3; (5) h = 0.5; (6) h = 1. The 
corresponding transition temperatures are: (1) T = 0.164; (2) T = 0.217; (3) T = 0.276; 
(4) T = 0.326; (5) T = 0.460; (6) T = 0.877. 

Fig. 3. Macromolecule order parameters: (a) w; (b) Si, (c) s 2 . Dependence on dimen- 
sionless temperature is shown for u = 0.8 and varying external field: (1) h = 0.01; (2) 
h = 0.1; (3) h = 0.2; (4) h = 0.3; (5) h = 0.5; (6) h = 0.8; (7) h = 1; (8) h = 2. The 
first-order transition temperatures are: (1) T = 0.143; (7) T = 0.715; (8) T = 1.764. For 
the values 0.01 < h < 1, the transition is a sharp crossover. 

Fig. 4. Macromolecule order parameters: (a) w; (b) si, (c) s 2 - Dependence on dimen- 
sionless temperature is shown for u = 1.5 and varying external field: (1) h = 0.01; (2) 
h = 0.5; (3) h = 1; (4) h = 2; (5) h = 2.5; (6) h = 3; (7) h = 5; (8) /i = 6. The transition 
temperatures are: (7) T = 4.547; (8) T = 6.4. For the values h < 5, the transition occurs 
as a sharp crossover. 

Fig. 5. First-order transition temperature T : (a) as a function of h at fixed u = 0.3; (b) 
as a function of w at fixed h = 0.5. 

Fig. 6. Qualitative illustration of the macromolecule swelling provoked by multiscale 
mesoscopic density fluctuations. 
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Figure 1: Macromolecule order parameters: (a) geometric weight of the strong-coupling 
phase w; (b) strong-coupling phase order parameter si, (c) weak-coupling phase order pa- 
rameter S2- Dependence on dimensionless temperature is shown for the repulsion parameter 
u = 0.3 and varying strength of external field: (1) h = 0.01; (2) h = 0.1; (3) h = 0.2; 
(4) h = 0.3; (5) h = 0.5; (6) h — 1. The corresponding first-order transition temperatures 
are: (1) T = 0.145; (2) T = 0.241; (3) T = 0.336; (4) T = 0.436; (5) T = 0.660; (6) 
T = 1.422. 



25 





Figure 2: Macromolecule order parameters: (a) w; (b) si; (c) s 2 - Dependence on dimen- 
sionless temperature is shown for the repulsion parameter u = 0.6 and varying strength of 
external field: (1) h = 0.01; (2) h = 0.1; (3) h = 0.2; (4) h = 0.3; (5) h = 0.5; (6) h = 1. The 
corresponding transition temperatures are: (1) T = 0.164; (2) T = 0.217; (3) T = 0.276; 
(4) T = 0.326; (5) T = 0.460; (6) T = 0.877. 
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Figure 3: Macromolecule order parameters: (a) to; (b) Si, (c) S2- Dependence on dimension- 
less temperature is shown for u = 0.8 and varying external field: (1) h = 0.01; (2) h = 0.1; 
(3) h = 0.2; (4) h = 0.3; (5) h = 0.5; (6) h = 0.8; (7) h = 1; (8) h = 2. The first-order 
transition temperatures are: (1) T = 0.143; (7) T = 0.715; (8) T = 1.764. For the values 
0.01 < h < 1, the transition is a sharp crossover. 
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Figure 4: Macromolecule order parameters: (a) w; (b) Si, (c) S2- Dependence on dimen- 
sionless temperature is shown for u = 1.5 and varying external field: (1) h = 0.01; (2) 
h = 0.5; (3) h = 1; (4) h = 2; (5) h = 2.5; (6) h = 3; (7) h = 5; (8) h = 6. The transition 
temperatures are: (7) T = 4.547; (8) T = 6.4. For the values h < 5, the transition occurs 
as a sharp crossover. 
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Figure 5: First-order transition temperature T : (a) as a function of h at fixed u = 0.3; (b) 
as a function of u at fixed h = 0.5. 




Figure 6: Qualitative illustration of the macromolecule swelling provoked by multiscale meso- 
scopic density fluctuations. 
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